NOTICE:  When  gpyernment  or  other  dravlngs,  s|>ecl> 
fieatlons  or  other  data  are  used  for  any  put^se 
other  than  in  conneGtlon  with  a  definiteiy  related 
goveroinent  proeurement  operation,  the  S« 
Government  therehy  incurs  n©  responslhllity,  nor  any 
obligation  '(diatsOever;  and  the  fact  that  '^e  Gksvern- 
ment  may  have  fOinmilated,  furnished.  Or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  penDission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  reiated 
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Eieveiosp'  new  radar  teefemqnes,  based  upon  stU'tistieali 
methods  j  whieh  are  applicable  to  automatic  radar  systemB 
with  inerease'd  detection  'capabilities  The  speeifie  phase 
reported  here  is  the  eonstruclion  and  statistical  evaluation 
of  an  electronic  generator  of  random  numbers  which  is 
used  in  laboratory  experiments,  evaluation  of  radar  detector 
performance  by  Monte  Carlo  methods,,  and  in  simulation  of 
physical  phenomena i 


RESULTS 


An  eleetronic  generator  of  random  numbers  has  been 
constructed  and  evaluated.  This  random  number  generator 
has  been  Used  to  generate  a  Sequence  of  ones  and  zeros 
With  a  predetermined  probahility  of  a  one  in  the  evaluation 
of  an  experi'mentai  sequential  detector.  The  device  has 
also  been  used  to  generate  random  sequences  of  ones  and 
zeros  in  simulating  the  input  to  an  experimental  track- 
before  ^detect  system  beihg  developed  at  NEL. 


Consider  the  RhfG  suitable  for  application  to  digital 
systems  requiring  the  use  of  random  numhers,. 
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Tfe  random  namfeer  generater  (ilMG)  desGritoed  in  this 
report  generates  14=^‘bit  random  namfeers  at  a  li^kc/s  ratei 
All  14  bits  aro  generated  siisinltaneonsly,  oaek  bit  feeing 
generated  ifey  a  separate  Sufe^^unit.  Tbe  il|MG  was  developed 
to  fee  nsed  in  eonjiunGtion  witfe  an  experinaental  seqnential 
deteitor  construGted  at  tfee  11*  Mavy  EleGtroniGs  Lafeora^ 
tory/^  A  seqnenee  of  ones  S-nd  zeros  witfe  a  prodetermined 
profeafeility  desGription  was  needed  as  an  input  to  tfee  de- 
teetor.  This  binomial  sequenee  was  generated  by  eompar^ 
ing  tfee  ©utput  of  the  fifefG  witfe  a  preset  I4^feit  binary  num^ 
ber^  iufeseqnently  tfee  BNG  output  has  been  used  m  the 
proeess  of  generating  random  sequences  of  ones  and  zeros 
in  simulating  tfee  input  to  an  experimental  traek^feefore- 
deteGt  System  being  developed  at  tfee  Navy  Electronics 
Laboratory.^  Tfeis  report  describes  the  generation  of  the 
random  numbers,  discusses  tfee  statistical  evaluation  of  tfee 
RNG,  briefly  describes  feow  the  first  application  mentioned 
above  was  made,  and  mentions  possible  otfeer  uses  of  such 
a  device  i 


Varioti^g  niaffiSricai!  Gatetalations  can  tee  carried  oat  wltte 
the  use  of  random  nambers.  The  term  '’Monte  Carlo''  is 
asaaliy  a;|)plied  te  any  sa eh  calculations .  As  an  exarhple. 


suppcise  it  is  d'esired  to  find  the  value  of  where 


6 

i  This  value  Can  tee  estimated  tey  setting  up  a 
stoChastiG  process  involving  a  fandom  vaf iatele  with  the 


property  that  S(y)  =  where  S{y)  is  the  expected 

O' 

value  of  l/i  Some  SOft  of  samplifig  procedure  Would  yield 
the  estimate.  One  method  of  sampling  would  tee  as  follows  j 
Ff  ohi  a  set  of  random  humbefS  uniformiy  diStribiUted  in  the 
interval  (0^  i)>  draw  two  sarnples  s  arid  t  ,  If  ^  =  1  for 

T_- 

/(s)  ^  ti  and  y  =  0  otherwise,  then  E'{y)  - 

b 


Monte  Carlo  methods  Can  also  tee  applied  to  the  evalu^ 
ation  of  e'haraCtef  istics  of  tests  of  statistical  hypotheses  . 
For  exarnple,  suppose  is  a  fandom  variateie  with  ffe=^ 
quenCy  functiOrii  /  {x,  9  ),  and  it  is  desifed  to  test  the  hypo^ 
thesis  Hq  :  f  ^  0o  against  the  alternative  :  6  =  0^  .  Ob- 
servations  are  to  tee  made  on  the  fandom  variable  x,  and 
some  criterion  is  eStabiished  by  which  either  ffo  Of  is 
aceepted  (e,  g. ,  the  Neyman*^Fearson  test).  Let  )  arid 

)tee  the  nurnber  of  times  out  of  m  tests  that  'Mq  and 
fespeGtively  are  accepted.  For  large  when  0  ??  0p,  then 
CL  (Mi  )  /m  I  and  when  '6  =  0i  then  ^  (Mq  )  /  ^ ,  whef  e  a 
and  ^  are  the  prob. abilities  of  type  I  and  type  II  errors,  re^ 
speGtively.  If  the  number,  n,  of  observations  required  in 
each  test  is  a  random  variable  (e.g.,  sequential  tests), 
and  if  N  is  the  total  number  of  observations  required  for  rn 
testa,  then  E  ( n  )< «:;  Nfm  . 

These  are  only  two  pi  many  applications  of  Monte 
Carlo  methods.  Other  applications  include:;  estimation  of 
collisions  involving  neutrons  or  other  particles,  solution 
of  differential  and  integral  equations,  and  the  inversion  of 
matrices , 


sEauENTiM  DETirraa  emuati0N 


The  seG^ndi  exaiHipie  m^htioned  afeove  suggested  evalu- 
Mifig  th^  eharaGteristics  of  lests  of  stafistiGal  hypotheses 
;by  Moiite  Carlo  rnethods .  It  was  for  this  purpose  that  the 
random  number  generator  was  developed^  &  sequential 
detector  to  detect  signal  in  the  presenee  of  noise  in:  rad^r 
video  was  to  be  evaluated.  The  video  was  quantized  into 
two  levels^  One  and  zeroj  hehce  the  deteetor  was  a  '’’bino¬ 
mial  sequeHtiai  deteGtori  "  ^  It  was  assumed;  that  when 
only  noise  waS  present  the  probability  of  observing  a  1  was 
and  when  :b©th  signal  and  riOise  were  preseht  the  proba¬ 
bility  of  observing  a  1  was  Px .  ^  Wald  sequehtial  prObabil^ 

ity  ratio  test^  was  to  be  Used  to  test  the  hypothesis  Mq  tp  -p© 
against  the  alterhate  :  p  -  Px  •  The  deteetor  eharaGteris  - 
tics  were  evaluated  by  Monte  Carlo  methods  utilizing  the 
random  humbejf  geherator.  This  required  gene  rating  a  se^ 
quence  where  X.  -  1  or  ^  .  sueh  that 

PrOfe  lx  •  1  i  -  n  where  p  is  the  parameter  of  the  bihOmial 

^  .....  ... 
distribution  GohSidered,  The  method  for  geheratihg  this 

Sequenee  will  be  described  later  i 

Briefly,  the  operatioh  Of  the  deteetor  was  as  follows: 

An  aGGumuiator  was  loaded  initially  with  a  quantity 
Dp  (O^Dq^I),  For  each  1  observed.^  a  quantity  D^{0^Dg^l) 
was  added  tp  the  contents  of  the  accumulator.  If  the  eOntentS 
of  the  aGcumuiator  were  theh  greater  than  or  equal  to  1, 
was  aecepted;  if  not,  another  observation  was  made.  For 
eaeh  0  observed,  a  quantity  was  subtracted 

from  the  contents  of  the  aeeumulator.  If  the  contents  of  the 
accunaulator  were  less  than  or  equal  to  zero.  He  was  aqeepted; 
if  not,  another  observation  was  made.  Hence,  after  initially 
loading  the  aeeumulator  with  Dq  ,,  observations  were  made 
until  one  oi  the  hypotheses  was  accepted.  By  repeating  this 
process  a  large  number  of  times  for  each  of  several  differ¬ 
ent  values  of  p  (p.  =  PrOb  ix^  -  l}:)  and  for  each  of  several 
different  values  of  Po,  P^ ,  and  Pj?,  the  detector  eharacter^ 
istics  were  evaluated,.  A  complete  description  of  the  detector 
and  its  evaluation  is  found  in  reference  2 .. 


A  'M(3©k  diagram  of  the  random  nnmfeer  generator 
(iRMiGrl  is  sfeowa  in  figure  1 .  The  Kil^G  generates:  a  14-bit 
binary  nnmfeer  with  the  property  that  the  probability  is  % 
that  any  partieular  bit  is  a  1.  Alse  each  bit  is  independent 
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Figure  i.  Block  diagram  of  the  rand  Cm  humber 
gene  rate r* 


of  ail  other  bits;  i.e,,  each  bit  is  generated  by  a  separate 
sub-unit.  The  output  of  the  RMG  is  either  seri&l  or  parallel 
depending  of  the  partieular  applieation  to  be  made  of  the 
output,  The  output  of  the  EKG  ean  be  eonsidered  to  be  a 
sequence  of  integers  uniformly  distributed  in  the  interval 
[Qj  2"^  -  1  ]  such  that,  for  each  observation  , 

Prob  for  h  =  ,0,  1,  . ,  ,  2^^  -1,  A  detailed  de- 

seription  of  the  ENG  follows. 

The  ENG  was  Gonstructed  from  commercial  plug-in 
type  digital  modules  together  with  laboratory  developed 
and  eonstructed  modular-type  pulse  generators.  The  digi¬ 
tal  modules  utilize  dynamic  circuitry  and  operate  synchro^ 
nously  with  respect  to  a  i-Mc/s  clock  frequency;  i.e.,  the 
"one”  state  is  a  l-Mc/s  train  of  -3  volt  pulses  and  the 
'  ■zero"  state  is  0  volt  dc .  The  modules  are  compatihle 


witl  inpHitS  from  other  modules  operatiBg  at  the  elock  fre- 
queiieyy  ahd  with  iflputs  Of  dc  voltages »  For  example^  a 
gate  with  a  do  iriput  of  -3  volts'  Will  assume  the  ’'©He''  state. 
Figure  2  shows  the  relationship  between  a  gate  output  and 
its  input.-  It  should  be  noted  that,  although  a  particular 
type  of  logic  device  is  discussed,  an  RNG  of  this  type  need 
not  be  restricted  to  their  usC  i 
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Details  of  Operation 

Figure  3  is  a  block  diagram  of  a  typical  subi-unit, 
fourteen  of  which  constitute  the  RNG.  Figure  4  is  a  dia^ 
gram  of  the  pulse  relationships  of  the  various  eomponenta 
of  a  typical  sub-unit,  figure  5  is  a  photograph  Of  a  thyra- 
tron  noise  generator  and  random^width-pulse  generator 
module,  and  figure  6  is  a  photograph  of  the  oomplete  RNG, 


PULSES  WITH  iNOlSE-MUDUlAlED  TRAILING  iDGES 


Figur»  $•  JBiook  diagram  of  a  typical  Muh~unit 


Figure  5*  Thy  rat  ran  noiae  generator  and  random- 
ioiid  th~  pul  a  e  g  en  e  ra  t  o  r  * 


Figure  6,  The  complete  random  number  generator. 
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The  following  desoription  refers  to  one  of  the  fourteen. 
Siiito  “Units  which  eonstitute  the  iKN'Gs  The  output  of  a  thyra- 
tron  noise  generator  is  arnplified  and’  USied  to  control  the 
pulse  width  of  a  pulse  generator »  The  output  of  the  noise  “ 
controlled“pulse  generator  is  a  pulse-width  naodulated  wawe 
of  positive  pulses  (referenced  to  “l®  volts)  with  their  lead¬ 
ing  edges  occurring  at  15  kc  /s.  The  pulse  width  varies 
randomly  feetween  fixed  limits,^  the  maximum:  width  depend¬ 
ing  on  the  maximum  modulation  applied.  The  pulse -width 
variation  about  the  Mean  is  approximateiy  Symmetric . 

The  pulse-width  moduiated  wave  is  applied  to  a  gate  which 
is  '''on"  during  the  inter  pulse  period  as  iadiGated  in  figure  4. 
Since  the  leading  edges  of  the  noiSe- modulated  pulses  occur 
at  fixed  intervals  and  the  pulse  widths  vary  randomly,  the 
iaterpulse  durations  also  vary  randomly,  ilence  the  output 
of  the  gate  to  which  the  pulse  "Width  modulated  Wave  is  ap¬ 
plied  is  a  succession  Of  random-length  packets  of  l-Mc/s 
pulses.  These  packets  trigger  a  hinary  flip-flop  Which 
Changes  state  foilowing  each  input  pulses  At  the  end  of 
each  packet  the  flip-flop  may  he  in  either  state,  one  or 
zero,  and  the  flip-flop  is  not  reset  hefOre  the  arrival  Of  the 
next  packet  of  pulses.  If  the  flip-flop  is  reset  between 
packets,  any  bias  toward  an  even  or  Odd  number  Of  pulses 
in  a  packet  shows  up  as  a  bias  in  the  number  of  zeros  Or 
ones  upon  sampling  the  flip-fiop.  In  order  to  smooth  out 
this  effect,  the  fLip-fiop  is  not  reset  between  packets-  it 
remains  in  its  prior  State  until  the  next  packet  arrives* 

An  analysis  of  this  smoothing  effect  is  given  in  the  appen¬ 
dix.  In  order  to  sample  the  state  of  the  flip-flop,  the  out^ 
put  of  the  flip-flop  is  AND- gated  with  a  synchronous  gen¬ 
erator  pulse  oecurring  at  15  kc/s  .  The  synchronous  gen¬ 
erator  pulse  is  delayed  in  Such  a  manner  that  it  appears 
before  the  flip-flop  stops.  This  procedure  is  used  to  avoid 
any  bias  which  might  result  from  a  tendency  of  the  flip-flop 
to  assume  a  particular  state,  partieularly  as  a  result  of 
"pulse  splitting"  which  might  oceur  on  the  last  pulse  of  a 
packet.  The  output  of  the  AND-gate  will  be  a  one  or  zero 
depending  on  the  state  of  the  flip-flop  at  the  time  of  arrival 
of  the  synchronous  generator  pulse ,  This  AND-gate  output 
is  the  output  of  one  sub-unit  of  the  KNGj  the  simultaneous 
output  from  all  14  sub-units  eonstltutes  the  RNG  output. 


Allhou-igii:  a  techuiqsae  for  generating  rando'ni'  binary 
digits  has  been  desoribe-dy  tbe  rantbod  ean  obviously  be 
extended  to  generate  nurnbers  to  any  integral  base  n  '(nag). 
The  randoin  width  pulses  are  applied  to  an  AbTD- gate  io- 
gether  with  the  output  of  a  clock  of  period  t  secondSy  the 
prodnct  tn  being  small  compared  with  the  smaliest  possible 
random  width  pulse.  The  output  of  this  gate  will  be  a 
series  of  paGfeetS  of  pulses,  each  packet  oontaining  a  ran¬ 
dom  number  of  clock  pulses  ..  Since  the  product  %t%  is 
small  with  respect  to  the  smallest  possible  random-widlh 
pulse,  the  probability  that  a  packet  of  pulses-  will  contain 
/d(mod  n)  pulses  (k  =  O;,  1,  . . .  yn  -  1)  is  approximately  l/n. 
Mence,  if  the  event  ''k  '(mod  n  )  pulses  in  a  packet is  called 
the  digit  h  (if-  0y  iy  . . n  ^  l,),  the  output  of  the  And -gate 
will  be  a  sequence  of  random  digits  with  approximately 
equal  probabilities  of  occurrence,  in  order  to  transform 
the  generated  sequenee  into  usabie  formy  the  AMD-gate  out¬ 
put  is  applied  to  a  counter  (mod  n )  and  the  state  of  the 
counter  is  sampied  in  the  interval  between  packets  of 
pulses.  The  n  possible  states  of  the  counter  will  be  identi¬ 
fied  with  the  digits  0,  ly  . . . ,  h  -  1  in  such  a  manner  that  0 
corresponds  to  0  (mod  n )  pulses  in  a  packety  1  corresponds 
to  1  (mod  h )  pulses  in  a  packety  etc .  Por  generating  ran^ 
dom  binary  digits  the  counter  takes  the  form  ©f  a  simple 
flip-flop  as  previously  describedi  for  decimal  digits  a 
decade  counter  would  suffice. 


Quantizing  the  RNG  Output 

As  mentioned  earlier,  the  motivation  for  developing 
the  ENG  was  to  evaluate  a  binomial  sequential  detector  by 
Monte  Carlo  methods,^  The  detector  input  was  required 
to  be  a  sequence  of  ones  and  zeros  with  a  predetermined 
(fixed)  probabiliiy  ©f  a  one.  This  detector  input  was  gen¬ 
erated  in  the  following  manner.  The  RNG  output,  A,  was 
compared  serially  in  a  comparator  (?,  with  a  preset  (fixed), 
number  Q'.  If  A  >  Q:,.  then  the  output  of  (7  was  a  one;  other¬ 
wise  the  output  was  a  zero,  Since  2^^  binary  numbers 


wer'e  eqpaily  likely  a>&  aft  output  fro-M  tke  RNG,,  the  fixe4' 
HU®i;feet  ‘Q  was  ehos'en  ifl  sueh  a  aiaanet  that  Proh  |i£*  -  1 } 
took  any  fixeh  value  wMeh  was  a  naiultiple  of  2“  and  Was 
hetween  zero  and  onei  The  output  of  the  oomparator  G  was 
taken  as  the  input  to  the  d'eteetof  * 


To  test  the  randomness  Of  the  ItNOi  four  tests  analog 
gous  to  those  used  hy  Kendall  and  Sniittf  were  applied  to  . 
the  output  of  eaeh  of  the  fourteen  suh^units,  Also  indepen= 
dence  hetween  sufe^units  wa$  tested  in  a  rnaaner  to  he  de- 
serihed  below. 

The  simplest  randoraikness  test  is  the  frequenGy  test* 
This  Gonsists  of  eounting  the  total  numiher  of  ones  in  a 
given  sanipiei  We  expeGt  the  nuniber  of  ones  to  he  One* 
half  of  the  total  sample  size. 

The  serial  test  ehecks  the  tendency  of  a  one  to  be  fol¬ 
lowed  by  a  one  or  by  a  zero,  etc.  This  test  is  performed 
by  sampling  in  pairs.  We  expeet  the  pairs  (1,  Ih  (Ij.  d)* 

(0,  l)i  and  (0,  0)  to  oeeur  equally  often. 

If  we  sample  the  output  of  a  sub-unit  in  blocks  of  three, 
for  example,  then  we  expect  one -eighth  of  the  total  number 
of  triples  to  contain  no  ones,  three -eighths  to  contain  a 
single  one,  three-eighths  to  contain  two  ones,  and  one- 
eighth  to  contain  three  ones  .  This  type  of  test  is  called 
the  poker  test. 


The  last  test  appiliedi  te  the  output  of  each  suh-uhit  is 
the  gap  tests  It  is  performed  by  ehoosing,  au  arbitrary 
bihary  digit,  say  a  one,  and  then  GOunting  the  number  of 
times  a  one  is  followed  by  a  one,  by  a  zero  and  then  a  one, 
by  two  zeros  and  then  a  one,  and  so  on.  In  other  words. 

We  eheck  zero  gap,  one  gap,  two  gap;  etcs  We  expect  zero 
gap  to  oGcur  one -half  the  time,  one  gap  to  oceur  one^fourth 
the  time,  etc  s 

Independence  between  sub -units  is  tested  by  s  imulta¬ 
neously  Sampling  the  Output  Of  two  or  more  of  the  sUb-unitS , 
If  the  outputs  of  two  sub-units  are  sampled  simuitaneously, 
then  we  expeGt  the  ordered  pairs  (0,  0),  (0,-1),  (1,  0),  and 
(1,  1)  all  to  OddUr  equally  Often.  This  is  sinailar  to  the 
serial  test,  except  that  here  each  member  of  the  observed 
pairs  is  the  ©utput  from  each  of  two  different  sub^unitSs  & 
similar  situation  results  if  three  or  more  of  the  sub“units. 
are  s ample d  simultaneously. 

A  further  test  of  independenGe  is  performed  by  quanti¬ 
zing  the  I4«bit  number  at  a  eertain  level  Q  as  previously 
deseribed.  Upon  choosing  Q,  the  probability  of  a  one  from 
the  output  of  the  quantizer  (comparator)  is  fixed  and  the 
expected  frequency  of  ones  is  known. 

The  Output  of  the  RN:G  was  also  indirectly  tested  for 
randomness.  As  mentioned  earlier,  the  quantized  output 
from  the  RNG  Was  used  as  the  input  to  an  experimental 
binomial  sequential  detector.  For  many  of  the  tests  per^ 
formed  on  the  experimental  detector,  exact  analytical 
oomputations  furnished  theoretical  results  ,  The  results  of 
the  Monte  Garlo  methods  were  in  excellent  agreement  with 
the  theoretical  results. 

The  test,  which  was  used  throughout,  tests  the 
goodness  of  fit  of  the  observed  frequencies  to  the  expected 
frequencies.  The  symbol  is  used  to  indicate  the  com¬ 
puted  x^  value  with  gj  degrees  of  freedom.  The  symbol  P 
is  used  for  Prob{  X^^Xg^  h  A*  is  used  to  indicate  the 
sample  mean.  The  statistic  for  the  serial  test  was  ehanged 
from  that  used  by  Kendall  and  Smith®  to  the  statistie  in¬ 
dicated  by  1,  J,  Good.  ® 
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RESULTS  IF  TESTS 

TaEbie  1  sbews  the  fesuits  of  the  serial  test  ©li  the  ont- 
put  ©f  ©fie  siah^^iUiMt  ©f  the  ran4©m  nuethef  generator.  The 
ex-pecte#  nnsnbeF  ©f  pairs  for  each  entry  in  the  table  is 
'25'00.,  The  data  for  the  fre^qnency  test  are  easily  'obtained 
frotn  the  sanie  datai  and  the  resnlls  of  this  test  are  in- 
eluded  in  the  table.  The  expected  number  ©f  ©bserved  ones 
in  each  row  is  5000,  The  Saniple  mean,,  p*,  is  given  for 
the  indiGated  GOlumns . 


Table  1 .  Serial  test  (pairs )  and  frequency  testi 
Ten  blocks,  l  O,  000  pairs  in  each  block. 


(1 ,1) 

Pa 

(1,0)' 

irs 

(0,1) 

(0'',Q> 

i 

Obs. 

:  Giles 

P 

2495 

2495 

2508 

0.040 

,84  ; 

4997 

e).oo4 

*.95 

2544 

2500 

:250G 

2440 

0 . 1  02 

.15 

:  9052 

1 .082 

.  30' 

2460 

2524 

2524 

2492 

0.92^ 

.34  ^ 

1  4934 

0,.  1 02 

•  79 

24:3 

2499 

2499 

2569 

0.002 

*96  ^ 

,  4932 

1,.'39o 

.17 

2502 

2502 

2457 

0.386 

.39 

^  5041 

0,.672 

*41 

24f8 

2539 

2539 

2424 

2.434 

.12 

9037 

0..  948 

.46 

2556 

2499 

2499 

2446 

0.002 

,96 

'  90  9  9 

1  ,210 

.27 

2512 

246,2 

246.2 

2964 

2,310 

.13 

!  4974 

0,270 

.60 

2547 

2465 

2469 

2923 

t  .960 

,  16  : 

;  9012 

0*093 

.81 

2467 

2496 

2496 

2941 

0 , 9.38 

,46 

4963 

0, 292 

.99 

24989 

24989 

24964 

0,049 

.<39 

90047' 

0.088 

.11 

,  2f0'58 

.  249.89 

, 24989 

, 24964 

. 90047 

Table  2  shows  the  results  of  the  poker  test  applied  to 
quadruples  from  the  output  of  one  sub-unit  of  the  random 
number  generator.  In  each  row  the  expected  number  of 
quadruples  containing  no  ones  or  four  ones  is  1024,  the  ex¬ 
pected  number  containing  a  single  one  or  three  ones  is 
4096,  and  the  expected  number  containing  two  ones  is  6144, 
Again  the  data  for  the  frequency  test  are  obtained  from  the 
same  data,  and  the  results  of  this  test  are  included  in  the 
table .  The  expected  number  ©f  observed  ones  in  each  row 
is  32,  768, 
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Table  2:*  Poiker  tes-t  (qBadrmples,)  and;  frequency  test » 
Ten  biocks>  16  ^  384  qnadfaiples  in  each  blocks 


total  10^88  4106’7  6ll8‘1  4t5:40  10264  1,807  -77  32760f  0.;O34  *85 

_p*  ,06279  *25065  *37342  .25049  *06265  .499886 

Table  3  shows  the  results  of  the  gup  test  appMed  to  the 
output  o!f  one  sub“un:it  of  the  randofn  numfeer  lenerater. 

The  expeeted  nuntber  ef  zero  gaps  is  8192,  the  expected 
number  of  one  gaps  is  4696;  etc.  No  frequency  test  results 
are  available  frona  these  data, 


Table  3,  Gap  test. 


Length  of 
'0 

gap 

1 

2 

3 

>4 

to 

P 

8238 

4102 

2046 

1006 

992 

1.585 

’h8l 

8175 

4118 

2097 

101'1 

983 

3. 132 

.55 

8169 

4134 

2035 

1036 

1  01 4 

0 .762 

.94 

8246 

4130 

2012; 

I1OI2 

984 

2.974 

*56 

8094 

'41  32 

2089 

1045 

1024 

2V740 

.60 

8041 

41 29 

2090 

1044 

1080 

7.364 

.12 

8187 

4091 

2048 

1037 

1 021 

0.174 

*99 

8117 

4083 

2060 

1r063 

1061 

3.620 

.46 

8194 

4072 

2065 

994 

1059 

2.357 

*67 

8205 

4057 

2045 

1020 

1057 

1.475 

..83 

Total  81662  41048  20587  10268  10275  1-757  *78' 

p*  .49843  .25054  .12565  .06267  .06271 
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Table  4  sbews  results  of  au  independence  test  between 
two  sub-units  of  the  randoiai  nunaber  generator  <  The  rnethod 
for  obtaining  the  data  is  described  previousiy.  The  ex¬ 
pected  number  of  pairs  for  each  entry  in  the  table  is  4096 . 

Table  4*  Independence  test  between  two  units » 

Ten  blocks,  16,  384  pairs  in  each  block. 


Pairs 


(  1:^1  ) 

(1^0)' 

(0^1) 

(;0>00  ' 

: 

41 23 

41 72 

3935 

41 04 

4.61 1 

,20 

4073: 

4162 

mo 

41 54 

4.707 

,20 

4106 

4151 

4021 

4106 

2.161 

.5< 

4068 

4054 

4073 

4189 

2.862 

,40 

41 56 

4025 

4063 

4140 

2,848 

,40 

4130 

4166 

40:24 

4064 

2,994 

• 

00 

4128 

4046 

4122 

4088 

1.041 

.80 

4199 

4041 

41 17 

4027 

4,598 

.20 

4031 

41 1  2 

4145 

4096 

1 ,680 

,65 

3983 

4092 

4250 

4059 

9,245 

,03 

Total  41002  41021  40790  41027  0.949  ,82 

p*  .2^026  .25037  ,24896  .25041 


Table  5  shows  results  of  quantizing  the  14-bit  binary 
numbers  from  the  random  number  generator.  The  quanti¬ 
zing  level  used  is  such  that  the  probability  of  a  one  is  the 
binary  fraction  .  01 01 01 0101 Q1 01  which  has  a  decimal 
equivalent  of  approximately  ,3333129883,  Hence  the  ex¬ 
pected  number  of  ones  in  each  block  is  5461  and  the  ex¬ 
pected  number  of  zeros  is  10^  923. 
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TaMe  5  .  ^antizatiori  ©f  14-feit  iiUBateer  , 

Ten  bloGka,  16,,  384  ©bseryations  in  each  feiocki 


Total 

P* 


Ohes  Zeros  P 


5366 

1  iors 

2.479 

*  12 

5412 

1  ©972 

©.659 

42 

534© 

1:1044 

4.021 

,..04 

5528 

108  56 

U213 

.27 

549© 

1©f94 

©.  231 

.63 

5446 

1  ©938 

©.  ©62 

,.  8© 

5547 

10837 

,2.031 

.15 

5472 

1  ©91 2 

©.©33 

..86 

5552 

10832 

:2.274 

.  13 

5466 

10918 

©.0069 

.93 

54619 

109221 

0, ©0222 

.96 

0333671  .6666321 


Table  6  shows  the  results  obtained  froni  Monte  Garlo 
eyaluation  of  the  binoHiial  sequential  detector  mentioned 
earlier,  for  the  particular  values  of  Dq  ,  D^,  and  Ey  in¬ 
dicated.  E(n)  is  the  expected  number  of  observations  for 
a  single  test,  L  ip)  is  the  probability  of  aceepting  the  null 
hypothesis  (Hq  )  when  p  has  a  particular  v^alue,  and  p  is  the 
probability  of  observing  a  one.  Data  for  the  columns 
headed  "Exact''  were  obtained  by  analytical  computations 
and  the  columns  labelled  "Exp,  "  contain  the  experimental 
results , 


Table  6 .  Ee)Suii:s-  of  Monte  Caf  lo  evaluation  of  the  binoHiial 
sequential  deteGtor  for  Dq  and  Og  -  l/S* 


p 

E(n.) 

Exact 

Exp, 

1 

Exact 

Exp.. 

.0312:5 

.  0&2:a 

18.0626 
23. 2574 

18.9150  I 
2:3.4211  ! 

1 . 00000 
:  . 99968 

1.00 00 
.  9997 

.  '00:375 

29. 8781 

30. 1512  : 

.99612 

,  99  5  7 

.125 

40.3776 

40. 1546  ! 

. 97390 

.  9744 

, 140625 

47.4225 

47.1290 

.94163 

.  9439 

. 15625 

55.2351 

55.2084  1 

1  .88108 

.  8787 

. 171875 

62.5432 

62 . 6890 

. 78143 

.  7819 

.1875 

67.310:2 

67.8328  1 

.  6  4:248 

.6483 

.203125 

67.8640 

68.0311  i 

i  .48309 

.  4729 

. 218  75 

64.1653 

64.2987  : 

:  .  33:333 

.  345  2 

.  234375 

57.8193 

56 . 9864 

.21499 

.  2204 

.  25 

50.6929 

50.9022 

.  132:64 

.  1319 

. 28125 

38.  2019 

38.2154  : 

;  . 04746 

.  0473 

.3125 

;  29. 4568 

29.4924 

!  .01662 

.0145 

. 34375 

23.5890 

23.7319 

1  .00583 

.  0064 

.  375 

19.5505 

19.5183  ' 

. 00205 

.0018 

. 40625 

16.6562 

16 . 8974  ' 

'  .00072 

.  0008 

The  results  of  the  statistical  tests  gave  no  reason  for 
rejecting  the  hypothesis  that  the  14-bit  numbers  generatod 
by  the  random  number  generator  are  uniformly  distributed. 
Also,  excellent  agreement  with  theory  was  obtained  in 
Monte  Carlo  evaluation  of  an  experimental  binomial  sequen-^ 
tial  deteetor- 
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AFPENUX;  mYSII  IF  ^*Nf)TiESETTIN@*^  PUP  ELOP 


We  assmtee  that  initially  the  flip^flep  is  in  the  zero 
state  ahd^  is  not  reset  after  being  triggered'  by  a  paGket» 

Let  p  be  the  iprobabllity  that  a  packet  has  an  odd  nnEnber  of 
puises  and  pj^  be  the  probability  that  the  ftip-flop  is  in  the 
one  state  after  having  been  triggered  by  k  packets  *  We  note 
that, p4  =  d  and p\=  pi 

Now  4 


The  general  solation  of  this  difference  eqaation  is 

p^  =  I  r  d(.l-:2p)^ 

Since  pg  -  0,  we  have  o  -  -§  » 


Therefore 


Let 


then 


and 


-  I  ( <  l-2p)^) 

p  -  I  ^  e,  |:>e  I  <  I 

p^  ^  |(l.(-2e)*) 


lini  p. 


_  1 
■'  2 


Thns  we  see  that  if  p  differs  frona  by  taking  k  large 
enough  P]^  will  be  as  close  to  ^  as  is  desired.  This  is  the 
motivation  for  not  resetting  the  flip-flop  after  each  packet. 
Closer  examination  of  the  above  analysis  reveals  however 
that  ''not  resetting"  is  not  a  "cure^nll,''  An  undesirable 
situation  exists  if  f  e  ,[  is  not  very  small.  If  we  assume  h 
large  enough  so  that  pj^^  then  Prob  { (L  1)  |  =  Prob  { (0,  0)} 
=  I  ( s  ^  ?)  and  Prob  {(1,  0) }  -  Prob  I  (0, 1)}  =  §(§  +  ?  ). 


TMs  M  f  e  f  IS  fiot  very  smatS  tile  .probabilities  as'S'ociated 
witb  pairs  will  be  severely  biased.  TO  take  an  extreme  ex*- 
ample  let  e  =  f*  Then  pjj.  =  §  (l  henee  =  0  and 

#2/^  +  1  ~  1-  Thus  the  pairs  '(1,  1)  and  (Ov  0)  can  never  ocGuri 

II  we  reset  the  fiip’-flop  after  the  arrival  of  eaGh  packet 
then  the  pair  probabilities  are : 

Prob  |  L)  I  -  i  +  e  ♦  e® 

Pr  ob  I  (1,  0  )'  I  =  ( I  +  e  .)  d  -  e  )  =  i  -  e® 

Prob  (0, 1)'  I  -  (I  -  e  )  (I  e )  =  I  >-  e® 

Prob  ,C  (0>  @1 1  -  d  -  ^  =  I  ^  e  ^  , 

©omparison  of  these  probabilities  with  those  of  the  "not 
reset"  situation  shows  that  two  of  the  pair  probabilities  are 
improved  by  resetting  while  two  have  been  biased  a  little 
morei 

Consideration  of  the  probabilities  associated  with 
n  ^tuples  shows  that  even  for  j:  e  f  small  some  of  the  proba'^ 
bilities  may  dev^iate  considerably  from  the  desired  values. 

To  obtain  some  idea  Of  the  bias,  a  long  count  of  Ones 
was  made  on  each  unit.  This  count  was  made  under  the 
reset  situation.  A  typical  result  for  one  unit  is 
j)*  -  .  499989140  based  on  a  total  saniple  of  5, 2^® .  The 
results  indicate  that  e  is  very  small. 


24 


INITIAL  DISTRIBUTION  LIST 


Chiefs  Bureau  of  Ships 
Code  320  Code  335  (3) 

GOd’e  360  ’(2)''CGde  6'7.0  {%) 

Code  684  {2^ 

Bureau  of  Naval'  Weapons’ 

DLi-3’  ©Li- 31 

RUDC-’2  RUDC-ll 

R-56 

Chief  of  Naval  Pfersonhel 
Library 

Chief  of  Naval  <^e rations' 

0tp-d7T  01^-94043 

0p-d3EG 

Chief  of  Naval  Research 
Code  455 

Commander  in  Chiefs  Pae  Fit 
Commander  in  Chief*  Lafit  Fit 
Commander  Operational  Test  & 
Evaluation  Force 
Deputy  Commander,  Operationai 
Te  st  ■&  Evaluafton  Fofcfe,  P'acif  iG 
Commandefi  Cruiser-Destroyer  Fori 
Pac  Fit  Lant  Fit 
Destroyer  Sevelopment  Groupi  Pac 
Commander  Training  Cornmand, 

Pac  Fit 

Commander  Amphibious  For, 

Pac  Fit  (2) 

Commander  Service  For,  Lant  Fit 
Naval  Air  ’Development  ’Center, 

Library 

Naval  Missile  Center 

Tech.  Library  Code  N32  32 
Naval  Air  Test  Center,  NANEP 
Naval  Ordnance  Lab.  i  Library 
Naval  Ordnance  Test  Station, 

Pasadena  Annex  Library 
Naval  Ordnance  Test  Station, 

China  Lake 

Naval  Weapons  Laboratory 
Naval  Radiological  Defense  Lab 
David  Taylor  Model  Basin 
Navy  Mine  Defense  Laboratory 
Navy  Underwater  Sound  Laboratory- 
Library  (3); 

ASW  Taeticai  School,  Lant  Fit 
Naval  Engineering  Experiment  Sta 
Library 

Naval  Research  Laboratory- 
Code  5330  Code  2027 
Cod.e  512P  Code  5400 
Naval  Ordnance  Laboratory,  Corona 
Navy  Underwater  Sound  Reference  Lab 
Library 

Air  Devei'ppment  Squadron  ONE  (VX-l) 
B.each  Jumper  Unit  ONE 
Beach  Jumper  Unit  TWO 
Office  ,pf  Naval  Research,  Pasadena 
Naval  Submarine  Base,  New  London 
Electronic  Officer 


Naval  Postgraduate  'School  - 
Library  (2 ) 

N  avy  Rep  r e  s  entative 
Project  LINGGL'N, 

Assistant  SECNAV 

Research  and  Development 
DODi  Reeearch  &  Engineering 
Technical  Library 
Army  Material  Command 
AMCRD'-DE-E-S 
-Assistant  'Chief  ef  StafL  G-2 
U  S  A^r  m  y,  'ID  B 

Aberdeen  Proving  Ground,  Tech  Lib,  (2l 
AM-STE-EL 

Army  Elect-rOnic  Proving  Ground, 
Technical  Library 

Redstone  ScientifiG  ilnfOrmation  Center 
Army  ElectroniGs  R&D  Laboratory 
SELRAySLr-ADT  SELR-A/SR 
Pieafihny  Arsefial 
Army  Research  Office  (Durham') 
Electronic  Developrnint  Activityi 
White  Sands 

Army  Transportation  Terminal 
Command,  Pac 

Army  Air  Defense  Board j  SUPPORT 
Deputy  Chief  of  Staff,  ©evelQp'ment,  USAF 
AFRST-ELycS 

Deputy-  Chief  o"f  Staffi  Operatipnsi  USA-F 
AFOCC-C/1 
Air  Defense  Command 
AP^AC-PL  AU0DA 
Air  Research  &  Development  Command 
SCSE 

Air  University,  Library  AyL3TT5028 
Alaskan  Air  Comrnand 
Electronics  (2)' 

Strategic  Air  Command 
OAST 

Air  Force  Missile  Test  Center 
MU^135, 

Rome  Air  Development,  Center 
RAALD 

Wright  Air  Development  Division 
ASAPRp.-pist  ASNVEG 
ASROQ-3 

Federal  Aviation  Agency 

NASA,  Langley  Research  Center  (3) 

Air  Weather  Service,  Scott  AFB 
National  Bureau  of  Standards,  Boulder 
Bureau  pf  Commercial  Fisheries 
Honolulu 

U,  S,  Weather  Bureau 
Washington  (2) 

Weather  Radar  Lab,  Oklahoma 
National  Security  Agency  (2) 

Georgia  Institute  of  Technology 
New  York  University 

Dept,  of  Meteorology  &  Oceanography 
University  of  Miami 

Marine  Laboratory  Library 


